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In the framework of the standard, hot big bang cosmological model the dynamics of the early evolution of the 
universe is controlled by the energy density of relativistic particles, among which neutrinos play an important 
role. In equilibrium, the energy density contributed by one flavor of relativistic neutrinos is 7/8 of that of 
the cosmic background radiation (CBR) photons. As the universe expands and cools, neutrinos decouple and 
their subsequent contribution to the energy density is modified by the relative heating of the CBR photons 
when electron-positron pairs annihilate. The small corrections to the post-e* annihilation energy density of 
the standard model neutrinos due to incomplete decoupling and finite-temperature QED effects are reviewed 
(correcting an error in the literature) and extended to account for possible additional relativistic degrees of 
freedom whose presence might modify the predictions of primordial nucleosynthesis and of the predicted CBR 
anisotropies. 



I. INTRODUCTION 

During the early evolution of the universe electroweak pro- 
cesses establish and maintain equilibrium between the stan- 
dard model neutrinos (e, fi, and r) and the cosmic background 
radiation (CBR) [y - In equilibrium the phase space dis- 
tributions are Fermi-Dirac for the neutrinos and Bose-Einstein 
for the CBR photons and the temperatures are equal (T 7 = 
T v = T Ve — T v — T Vt ). The energy densities are related by 



Pv = Pv e 



/V = zPi- 



(1) 



During early, radiation dominated (RD) epochs the universal 
expansion rate is uniquely determined by the energy density in 
relativistic particles. Prior to e ± annihilation, at a temperature 
of a few MeV, the energy density receives its contributions 
from CBR photons, e pairs, and three flavors of neutrinos, 
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P = P~,+ Pe+iPv = -£"P 7 . 



extra energy density by normalizing it to that of an equivalent 
neutrino flavor [01, 



p x = AN„p„ = 7 -AN uPl . 



(4) 



In the following it is assumed that X has decoupled prior to 
e^- annihilation and its temperature (or, equivalently, its num- 
ber density) is used to define a "comoving volume". For each 
such "neutrino-like" particle {i.e., a two-component fermion), 
if T x = T„, then AN U = l ; if X is a scalar, AN„ = 4/7. 
However, it may well be that X has decoupled earlier in the 
evolution of the universe and has failed to profit from the heat- 
ing when various particle-antiparticle pairs annihilated (or un- 
stable particles decayed). In this case, the contribution to 
AN V from each such particle will be < 1 (< 4/7). Since we 
are interested in the deviations of T 7 and T v from Tx when 



e ± pairs annihilate, for convenience and without loss of gen- 
erality, we may define Tx = T v prior to e ± annihilation. In 
the presence of this extra component, the pre-e* annihilation 



(2) energy density in eq. || is modified to, 



At this time (T ~ few MeV) the neutrinos are beginning to de- 
couple from the photon - e ± plasma and the neutron to proton 
ratio, which is key to the primordial abundance of 4 He, is de- 
creasing. As a result, the predictions of primordial nucleosyn- 
thesis depend sensitively on the early expansion rate. During 
these RD epochs the age and the energy density are related by 
^-Gpt 2 = 1, so that the age of the universe is known once 
the temperature is specified, 
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(5) 



The extra energy density speeds up the expansion of the uni- 
verse so that the right hand side of the time-temperature re- 
lation in eq. || is smaller by the square root of the factor in 
parentheses in eq. || 



t T 7 = 0.738 MeV 2 s. 



(3) 



II. COMPLETELY DECOUPLED NEUTRINOS 



In many extensions of the standard models of cosmology 
and of particle physics there can be "extra" energy density 
contained in new particles or fields, px- When X behaves 
like radiation, that is when its pressure and energy density are 
related by px = \px, it is convenient to account for this 



For standard model neutrinos the electroweak interaction 
rates drop below the universal expansion rate when the uni- 
verse is less than one second old and the temperature is a 
few MeV. The "standard", zeroth-order approximation that the 
neutrinos are completely decoupled prior to annihilation is 
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sufficiently accurate for most cosmological applications. In 
this approximation the CBR photons get the full benefit of the 
energy/entropy from the annihilating e^- pairs, maximizing 
their heating relative to the decoupled neutrinos (and to any 
other decoupled particles, such as X). After e^- annihilation, 



T v = T X = ( — ) 1/3 T 1 = 0.7138T 7 



(6) 



In this approximation the Fermi-Dirac phase space distribu- 
tions of the decoupled neutrinos is preserved as the universe 
expands and p„ e = p Vft = p Ur = p v (and px = AN v p v ), 
where 



pv_ 

Pi 



1(1jL\i - Z(_lW3 = n 2271 



and 
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(7) 



(8) 



In this zeroth-order, full decoupling approximation, the post- 
e ± annihilation relativistic (i?) energy density is, 

— = 1 + 1(tt) 4/3n » = 1-6813(1 + 0.135lAAy, (9) 
p 7 oil 

where N v = 3 + AA^. As long as the universe remains radia- 
tion dominated, the age and the photon temperature are simply 
related by, 



t°T* = 1.32(1 + 0.1351AA,,)- 1 / 2 MeV 2 s. 



(10) 



Both pR and the time-temperature relation play important 
roles in establishing the amplitudes and angular scales of the 
CBR anisotropies and any corrections to these zeroth-order 
results may bias the interpretation of the precise data from 
current ground-based and future space-based telescopes. 



III. PARTIALLY COUPLED NEUTRINOS 

It has long been known that the neutrinos are not entirely 
decoupled from the electron-positron-photon plasma during 
e ± annihilation [Q], As a result, the neutrinos get to share 
some of the annihilation energy with the photons and the 
relative heating of the photons is reduced from the zeroth- 
order approximation estimate. Furthermore, since the weak 
interactions coupling the neutrinos and the electrons are en- 
ergy dependent, the neutrino phase space distributions are dis- 
torted by preferential heating of the higher energy /momentum 
neutrinos and the resulting distributions are no longer Fermi- 
Dirac. In a series of increasingly detailed calculations, many 
authors have tracked this evolution [^] - JT2| | . Although there 
are small differences in the quantitative results which may be 
traced to the differing approximations, the overall agreement 
among them is excellent. The results of the detailed and ex- 
tensive calculations of Gnedin & Gnedin [O], which cover 



some seven orders of magnitude in the neutrino momentum, 
are adopted here. 

Recall from eq. 7 that in the fully decoupled approxima- 
tion, Pv/p-y — 0.2271. Since the electron neutrinos participate 
in charged-current as well as in neutral-current weak interac- 
tions, they remain coupled longer than do the /jorr neutrinos 
and, thereby, are heated more. In following the decoupling 
carefully, Gnedin & Gnedin find, 



^ = 0.2293 = -(— ) 4 / 3 x 1.0097, 
p 1 8 11 



while 



Pi 



^ = 0.2285 = l(— ) 4 / 3 x 1.0061. 
Pj 8 11 



(11) 



(12) 



Further, because the energy/entropy from the e is now being 
shared by the neutrinos, the heating of the photons relative to 
the X (assumed to be fully decoupled prior to annihilation) 
is reduced, so that Tx = 0.7144T 7 (compared to the zeroth- 
order estimate in eq. 6) and, 



7 ,T 



7 , 4 



Px_ = 

Pl 8 v T 7 ' ' S ' 1 1 



= ) 4 AA„ = -( — ) 4 / 3 AA„ x 1.0036. (13) 



We are now in a position to combine the post- 
e ± annihilation energy densities of the photons along with 
that of the neutrinos (and possible As) to find the relativistic 
energy density in the incompletely decoupled (ID) neutrino 
approximation. This result may be written in analogy with 
eq. 9 by replacing N v with N^ D , 



n ID 7 4 

p 7 8 11 



where 



N, 



N, 



ID 



ID 



3.022 + 1.0036AAT,,, 



N v + (0.022 + 0.0036AAy 



(14) 



(15) 



(16) 



Thus, for the standard model case of N u = 3, the post- 
annihilation energy density corresponds to 3.02 "equiva- 
lent" neutrinos. Although all authors ( [|J - [[j~2]]) agree quanti- 
tatively with this AiV„ = correction, a different value (3.03) 
appears in the Lopez et al. paper [|l3|] . It is difficult to iden- 
tify the source of the Lopez et al. correction, 8Nl° = 0.03, 
which differs from that found earlier by some of the same au- 
thors [|8 14]. The Lopez et al. value, modified by the QED 
effect to be addressed next (SN^ ED = 0.01), seems to have 
propagated in the recent literature (see, e.g., \ 15,1^]) and, in- 
deed, 3.04 is the recommended default value for N u in the 
CMBFAST code of Seljak & Zaldarriaga [0, widely utilized 
in analyses of the CBR fluctuation spectra. The second term 
on the right hand side of eq. |l6[ the difference between Nl D 
and N v , is the incomplete decoupling correction to the stan- 
dard model result, including the effect of extra relativistic en- 
ergy. However, there is another effect which, while small, is 
the same order of magnitude as this correction for incomplete 
decoupling. 
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IV. QED CORRECTION 



In the zeroth-order approximation, the 7 



p 7 + 8 l ll J v ' 



7 , 4 



(21) 



: plasma is 

treated as a gas of free, non-interacting particles prior to 



where eq. ^Ojprovides the connection between N*?f and N v . 
For the standard model case of N v = 3, the correction to the 
e ± annihilation When finite temperature QED corrections post-e* annihilation energy density, 0.034, is small; but, for 

AA„ > 0, this correction term grows. The new term, propor- 
tional to AN V , will dominate the correction for AA„ > 4.5, 
and should not be ignored if precision requires that the first 
term be included. The difference between the zeroth-order 
energy density (eq. 9) and the corrected one (eq. 21) is also 
shown in Figure 1. It is interesting to note that in a joint 
BBN + CBR analysis of the constraints on N?*, Hansen et 
al. ( [[20]]) recommend a standard model (AA„ = 0) value of 
iV^" = 3.034, in excellent agreement with that presented 
here, but they make no mention of the additional corrrection 
when AN,, differs from zero. 



are included ( [ 18 19|]), the energy density and pressure of the 
7 - e ± plasma are reduced. As a result, when pairs an- 
nihilate they actually have less entropy to share with the pho- 
tons (and the incompletely decoupled neutrinos), than would 
be estimated neglecting this QED correction. Accounting for 
this reduction leads to fewer CBR photons (and neutrinos) in 
the post-e ± annihilation comoving volume, corresponding to 
a lower photon temperature. To quantify this effect ( Jl8| , |l9| l) 
it is convenient to compare the post-e 1 * 1 annihilation photon 
temperature to that of the fully decoupled X particles. 



(Tx/T^qed = (4/ll) 1 / 3 (l + 9.6 x 10- 4 ). 



(17) 



Thus, in comparing the total energy density in relativistic par- 
ticles to that in CBR photons alone, the relative contribution 
of the neutrinos (and Xs) is enhanced with respect to the con- 
tributions in the zeroth-order approximation (eq. 9) and in the 
incompletely decoupled approximation (eq. 14). To account 
for the finite temperature QED correction, the (4/11) 4 / 3 fac- 
tors on the right hand sides of eqs. 9 and 14 should be re- 
placed by (Tx/T 1 )q ED , leading to a new, "effective" number 
of equivalent neutrinos N^ * , 

Nt ff = (ll/4) 4 / 3 (T x /T 7 ) 4 D Nl D = 1.Q03&/V™, (18) 



so that 



or 



N, 



eff - 



3.034+ 1.0074AA„, 



NfJ f - N v = 0.034 + 0.0074AiV„. 



(19) 



(20) 



This total correction to the standard model, zeroth-order ap- 
proximation, N v = 3 + AN U , is shown in Figure 1 as a func- 
tion of any extra energy density, measured by AN V . 



V. DISCUSSION 

Prior to e^- annihilation the total energy density in the stan- 
dard model (N v = 3) is given in terms of the CBR energy 
density by eq. 2. In the presence of extra (relativistic) energy 
density, AA„ > 0, the pn — p 1 relation is modified to that 
in eq. 5. After e annihilation the relativistic energy den- 
sity is shared by photons and decoupled neutrinos (and As) 
and the total energy density is modified by the relative heat- 
ing of the photons, neutrinos, and As. When careful account 
is taken of the incomplete decoupling of the neutrinos during 
annihilation, as well as of the finite temperature QED ef- 



fects on the equation of state of the pre-annihilation plasma, 
the total energy density (p^ 
photon energy density, is 



= p R ), normalized to the CBR 




FIG. 1. The solid curve is the correction to the zeroth-order result 
for the equivalent number of neutrinos as a function of the "extra" 
equivalent number of neutrinos. The dotted curve shows the frac- 
tional change in the post-e* annihilation age when the universe is 
radiation dominated. The dashed curve is the difference between the 
corrected, and zeroth-order radiation energy densities in units of the 
CBR photon energy density. 

Since the post-e ± annihilation expansion rate is determined 
by pn (as long as the universe remains RD), the age - tem- 
perature relation will be modified from its zeroth-order form 
(eq. 10), 



t° 1 

7 = ' 



+ 0.1351(A^-3) 
1 + 0.1351(^-3) 



1/2 



(22) 



The fractional change in the age of the universe (At/t = 
(t — t°)/t) is shown in Figure 1. Since the CBR fluctuation 
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spectrum is being painted on the microwave sky during those 
epochs when the universe is making the transition from radi- 
ation to matter dominated, modifications to the zeroth-order 
expressions for the expansion rate and the radiation energy 
density will effect the details of the resulting anisotropics. 

As cosmology enters this new era of precision science, it 
may be necessary to account for even the very small changes 
from the zeroth-order approximation summarized here. For 
example, Lopez et al. ( [13]]) suggest that measurements ac- 
curate to SN V » 0.03, a ~ 1% determination, will be pos- 
sible. At present, such precision seems a distant dream. For 
example, for AiV^ = 0, the effect on the BBN-predicted he- 
lium abundance is very small ( [ p"4[|I^ ]), 5Yp = 1.5 xl0~ 4 , a 
correction buried in the overall uncertainty in the BBN pre- 
diction (X, 4 x 10~~ 4 [21]). The small speed up in the post- 



between Ny 1 and N v at the 1% level are beyond our grasp. 
We can only hope that our experimental colleagues will rise to 
the challenge and provide data of such exquisite precision that 
the corrections reviewed and summarized here will be impor- 
tant. 
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